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I .  INTRODUCTION 


Let  X| ,  X2>  . ...  Xn  be  independent,  identically-distributed  random 

variables.  Then  Yj^  <  Y2  <  ...  <  Yn>  where  the  Y^ * s  are  the  X^s  rearranged 

in  order  of  increasing  magnitudes,  are  defined  to  be  the  order  statistics 
corresponding  to  the  original  random  sample.  Order  statistics  find 
immediate  application  in  the  design  and  evaluation  of  logical  structures 
which  make  decisions  based  on  the  relative  values  assume-;  1/  a  set  of  n 
random  variables.  One  particularly  interesting  app-ii-.'  ion  involves  the 
determination  of  time  windows  for  firing  impulses  in  a  fuzing  system. 

This  is  the  problem  which  motivated  the  work  on  which  we  are  reporting. 
However,  the  work  is  general  in  nature  and  may  prove  useful  for  other 
applications  in  the  armament  research  and  development  community. 


II.  STATEMENT  OF  THE  PROBLEM 

In  the  particular  problem  which  we  addressed,  a  fuze  contains  N  detonators, 

K  of  which  must  function  within  a  specific  time  span.  Furthermore,  the  second 
detonator  (which  functions  at  time  Y2)  partitions  the  time  span  into  two  sub¬ 
intervals.  The  first  subinterval  [Y2-  5j,  Y2]  is  examined  to  determine  if  the 

first  detonator  functioned  within  that  time  segment,  and  the  second  subinterval 
[Y2,  Y2  +  <$2]  is  monitored  to  count  the  number  of  additional  detonators  activated 

during  that  period  of  time.  If,  within  the  time  interval  [Y2  -  6j,  Y2  +  <$2], 

K  detonators  have  functioned,  then  the  command  to  fire  will  be  initiated; 
otherwise,  it  will  not.  The  times  to  function  for  the  detonators  are  random 
variables  and,  as  such,  can  be  characterized  by  a  cumulative  distribution 
function  F.  Assuming  that  the  time  to  function  of  each  detonator  is  identi¬ 
cally  distributed,  then  the  problem  consists  of  expressing  the  probability 
of  fuzing  as  a  function  of  K,  N,  fi2,  and  p* 


III.  SOLUTION 


Let  be  the  time  to  function  of  detonator  i  in  its  operating  environ¬ 
ment.  Then  X^,  X2,  ...,  X^  are  independent,  idei  tically-distributed  random 
variables;  and  we  can  define  Y^,  Y2,  ...,  Y^  to  je  the  order  statistics 
corresponding  to  the  X^'s,  For  our  problem,  if  Y2  -  Yj  <  6^,  we  are  interested 
in  the  probability  that  Yj,  -  Y2  *  ^2*  however,  ^  Y2  "  Y1  >  ^1’  we  need  to 
determine  the  probability  that  YR+1  -  Y2  <  62  assuming  K  +  1  <  N.  Therefore, 
we  need  to  evaluate 


Pr  (warhead  fuzing}  ■  Pr  (Y2  -  Yj  <  5j}  Pr  {Y^  -  Y2  <  <S2  \  Y2  -  Y1  <  5j} 


+  Pr  (Y2  -  Yj  >  fij}  Pr  (YK+1  -  Y2  <  62  |  Y2  -  Yj  >  «jJ  . 


(1) 


7 


which  is  equal  to 


Ji  ■  (h-k) r i.  fM  f!1  f(v*w)  1F(U*IV. V*“) '  ’’(''♦"’i"'5 

•  f(u+v+w)  [1  -  F(u+v+w)]^"^  du  dv  dw  .  (8) 

In  an  analogous  manner  we  can  obtain  the  joint  probability  density 
function  of  the  1st,  2nd,  and  K+lst  order  statistics;  and,  letting 
u  ■  y^+2  ~  V 2*  v  "  ^2  "  ^1  an^  w  ■  we  can  obtain  the  necessary 

probability  for  the  second  term  on  the  right-hand  side  of  Equation  3. 

That  probability  is  equal  to 

4*00  +00  j 

J2  ■  TTSTtfepT  J_  «M  |  f(v«)  J,2[ Ptu*v«o  -  P(v«0]K-2 

•  f(u+v+w)  [1  -  F(u+v+w)]N"K_1  du  dv  dw.  (9) 


The  probability  of  fuzing  is  then  just  the  sum  of  Equation  8  and  Equation  9. 

Appendix  B  contains  a  computer  program  which  evaluates  these  inte¬ 
grals.  In  its  current  form  it  assumes  that  the  distribution  of  the 
times  to  function  of  the  detonators  is  normal;  however,  it  can  be  easily 
modified  to  change  this  distribution.  The  program  requires  as  input  N, 

K,  dj,  fi2,  and  a  (standard  deviation  of  the  assumed  distribution). 

IV.  RBSULTS 

For  the  problem  we  addressed,  the  value  for  a  was  specified  to  be 

10“S  seconds.  Figure  1  presents  the  results  of  changing  and  52  for  a 

fuzing  system  in  which  K  is  equal  to  N-l.  In  Figure  2  we  considered  a 

fuze  with  eight  detonators;  and,  keeping  a  equal  to  10-5,  we  varied  K  as 
well  as  6^  and  62<  Of  course,  given  any  four  of  the  five  input  variables 

(o,  K,  N,  Sj,  and  6^),  we  can  obtain  a  specific  probability  of  warhead 

fuzing  by  parametrically  varying  the  remaining  variable. 


V.  APPLICATION  TO  FUZB  DESIGN 

An  important  consideration  is  to  determine  optimal  and  62  based  on 

a  specified  cj,  K,  and  N  and  a  desired  probability  of  warhead  fuzing.  To  do 
this  several  factors  must  be  considered.  To  maintain  the  reliability  of  the 
fuze  it  is  necessary  that  the  interval  [Y2  -  6j,  Y2  +  62]  be  sufficiently 
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FIGURE  1.  PROBABILITY  OF  WARHEAD  FUZING  IN  AN  N-l  OUT  OF  N  FUZING  SYSTEM  (o  * 


FIGURE  2.  PROBABILITY  OF  WARHEAD  FUZING  IN  A  K  OUT  OF  8  FUZING 


large  to  allow  for  the  occurrence  of  K  detonations.  However,  safety  must 

also  be  considered;  and  here  it  is  desirable  that  the  interval 

[Y2  -  Y 2  +  djl  be  as  small  as  possible  to  preclude  the  occurrence  of 

spurious  detonations.  These  conflicting  goals  could  be  accomplished  through 
a  solution  to  the  following  optimization  problem: 

Given  a,  K,  and  N,  determine  S2  >  0  for  which 
+  62  minimum  subject  to  Jj  +  Jj  >  a  (reference 

Bquation  8  and  Equation  9)  where  a  is  the  desired 
probability  of  warhead  fuzing. 

Consideration  of  this  problem  is  beyond  the  scope  of  this  report.  Here, 
we  will  content  ourselves  with  the  determination  of  reasonable  initial  values 
of  6 j  and  <S2  £°r  use  in  the  computer  algorithm.  Such  values  may  then  be 

adjusted  as  required  to  obtain  the  desired  probability  of  warhead  fuzing. 

For  this  purpose  practical  initial  values  for  6^  and  S2  would  be  B[Y2  -  Yj] 

and  B[Yk  -  Y2]  respectively,  where  E[.]  is  the  expected  value  of  a  random 

variable. 

For  the  general  case  we  can  define  WRg  ■  Yg  -  YR,  and  from  consideration 
of  the  material  in  Appendix  A  we  obtain  the  probability  density  function 

00 

^wrs^  "  fcys  *  *  ^rs  ^ 

•  [F(x  ♦  wRg)  -  F (x) ] s~^“ 1  f(x  +  wRS) 

•  [1  -  F(x  ♦  wrs)]N"S  dx  (10) 

where  CRS  -  (r.i) |  (S-R-l) I  (N-STT  1 
To  determine  E[WRg]  ■  E[Yg  -  YR]  we  must  evaluate 

4.  CO 

BfWRS^  “  ^0  WRS  f{v'RS-)  dwRS *  (11) 

Therefore,  we  can  obtain  E[Y2  -  Yj]  by  evaluating 

+00  +00 

E[Wi2]  -  C12  /8  /  w12  f(x)  f(x  +  w12) 

•  [1  -  P(x  +  w12)]N“2  dx  dw12;  (12) 
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and  in  an  analogous  manner  we  can  obtain  E [YK  -  Y2]  by  evaluating 


E[W2](]  ■  C2|(  J  ^  w2K  F(x)  f(x) 

•  [F(x  +  w2K)  -  F(x)]K"3  f(x  +  w2K) 

•  [1  -  F(x  +  w2k5]N"K  ^x  dw2K  ’  O3) 

The  pth  quasi-range  of  N  sample  values  is  defined  to  be  the  range  of 
<M  -  2p)  values,  omitting  the  p  largest  and  the  p  smallest.  For  the  special 
case  in  which  K  ■  N-l,  W2K  ■  yn_1  -  Y2  becomes  the  first  quasi-range;  and 

this  distribution  has  been  investigated  and,  to  some  extent,  tabulated. 

For  example,  H.  Leon  Harter*  has  provided  the  expected  values  of  the  first 
quasi-range  »  E [W 2  N  j]  for  samples  from  a  normal  distribution  with 

mean  zero  and  variance  one  and  for  values  of  N  from  four  to  one  hundred. 
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APPENDIX  A 

DERIVATION  OF  THB  JOINT  PROBABILITY  DENSITY  FUNCTION 


As  shown  by  H.  A.  David,  we  can  obtain  the  joint  probability  density 
function  for  Y1(  Y2,  and  YR.  Let  XJ(  X2,  Xn  be  a  random  sample  from 

a  population  with  density  f(x)  ■  F'(x)}  and  let  Yj,  Y2,  ...»  Yn  be  the 

corresponding  order  statistics.  Then  the  probability  density  function  for 
the  rth  order  statistic  may  be  derived  by  considering  the  following  config¬ 
uration: 


That  is,  <  x  for  r-1  of  the  X^,  x  <  X^  <  x  +  Ax  for  one  X^,  and 
X^  >  x  +  Ax  for  the  remaining  n-r  of  the  X^.  The  number  of  ways  this 
combination  of  events  can  occur  is 


and  each  such  way  has  probability 

[F(x)]r-1  [F(x+Ax)  -  F(x) ] 1  [1  -  F(x+Ax) ]n”r. 


Therefore,  we  have 

Pr  {x  <  Yr  <  x  ♦  Ax}  -  F:T)f  (n.r)r 
•  [F(x) ]r-1  [F(x+Ax)  -  F(x) ]  [1  -  F(x+Ax)]n"r  +  0(Ax2)  (Al) 

2  2 

where  0(Ax  )  means  terms  of  order  (Ax)  and  includes  the  probability 

of  realizations  of  x  <  Y„  <  x  +  Ax  in  which  more  than  one  X.  is  in 

r  x 

(x,  x  +  Ax) .  Dividing  both  sides  of  Equation  A1  by  Ax  and  then  letting 

Ax  •+■  0,  we  obtain 


fY  (x) 

r 


[FCx)]r_1  f (x) [1-F(x)]n"r. 


David,  H.  A.,  Order  Statistics,  John  Wiley  and  Sons,  New  York,  1970,  Pp.  9-10. 


In  a  similar  manner  we  can  derive  the  joint  probability  density 

function  of  Yr  and  Ys; 


r-l 

s-r-1 

n-s 

X 

x+Ax  y 

y+Ay 

and  through  an  analogous  argument  we  obtain 

V  '*•>”  ■  TnTrrrWrtn^T  nwiMfw 

X  s 

•  [F(y)  -  Fto]3"1’1  f(y)  [1  -  F(y)]n'S.  CA3) 

Finally,  for  the  joint  probability  density  function  of  Yr,  Yg,  and  Yt 


r-l 

1 

s-r-1 

1 

t-s-1 

l 

n-t 

X 

x+Ax  y 

y+Ay  i 

i+Az 

and 

fYr,Ys,YtCx,y,z)  "  ‘(r-l)  I  (s-r-1)  1  (t-s-1)  I  (n-t)  I  [F(x)]r 
•  [F(y)-F(x)]S"r~1  f (y)  tFCzJ-FCy)]1"9’1  f(z)  [l-PCt)]n_t.  (A4) 

For  the  case  »«1,  s»2,  and  t-k,  we  obtain 

fYi,Y2,Yk(x*y'z)  "  -nrsTrin-kyi fCx)  fCy) 

•  [F{z)-F(y)]k“3  f(t)  [1-F(0]n'k.  (A5) 
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APPENDIX  B 
COMPUTER  PROGRAM 

We  are  presenting  here  the  computer  program  which  evaluates  the  desired 
probability  function.  As  noted  in  the  program  comments,  we  have  assumed  that 
the  times  to  function  of  the  detonators  are  normally  distributed  with  mean 
2 

zero  and  variance  a  .  However,  with  just  a  few  changes  to  the  subroutines, 
a  different  distribution  may  be  assumed.  To  do  this,  all  cards  containing 
'FORMAL"  in  columns  74  through  79  must  be  replaced  by  others  with  the  approp¬ 
riate  probability  density  function  or  cumulative  distribution  function. 


UUUUUUUUUUUVUUUUUUUUUUUUU  O  O  o  u  u  oooo 


PROOPAM  OC TON  <INPWT.0UTMUT.TAPES»INPUT,TAPE6»0UTPUT> 


THIS  PROORAM  DETERMINES  THE  PROBABILITY  THAT  A  OUT  OP  N 
OETONATORB  PILL  FUNCTION  WITHIN  A  OIVEN  T IMESPAN.  K  MUST 
•I  ORCATCR  THAN  OF  EQUAL  TO  It  ANO  K  MUST  BE  LESS  THAN  OR 
EQUAL  TO  N*  THE  OISTRIRUTION  OF  THEIR  FUNCTIOnInO  IS  ASSUMED 
TO  BE  OAUSSIAN  KITH  A  MEAN  EQUAL  TO  OtO  AND  A  STANDARD 
DEVIATION  EQUAL  TO  SIONA,  THE  SOLUTION  IS  DERIVED  THROUQH 
THE  USE  OF  ORDER  STATISTICS  AND  IS  OBTAINED  BY  EVALUAT INO 
A  TRIFLE  INTERNAL. 

IF  THE  NUMBER  OF  OETONATORB  IS  BREATER  THAN  THIRTY. 

THEN  THE  TOLERANCE  LIMITS  OF  THE  INTEGRALS  MUST  BE  REDCFINEO. 
THAT  IS.  THE  VARIABLE  ’ERROR'  IN  THE  MAIN  ROUTINE  AND  THE 
VARIABLES  'TOLXtTOLV '  IN  SUBROUTINE  'DIST'  SHOULO  BE  AOUUSTED. 
with  the  TOLERANCE  LIMITS  CURRENTLY  IN  THE  PROORAM. 

THE  RESULTINO  FROSARILITIES  ARE  CORRECT  TO  TWO  DECIMAL  FLACES. 

INPUT  IS  AS  FOLLOWS  .... 


N  . . NUMBER  OF  DETONATORS 

K  .................  NUMBER  OF  DETONATORS  TO  FUNCTION 

SI OH A  .............  STANDARD  DEVIATION  OF  DISTRIBUTION 

OELTI  .  TIMESFAN  (FIRST  TO  SCCONO  OETONATORB) 

OElTS  . . .  TIMESFAN  (SECOND  TO  LAST  DETONATORS) 


REAL  NFACT .K2FACT.K3FACT .NKFACT.NKiPACT 
DIMENSION  ICRRIB) 

COMMON  /COM1/  N.K.SIBMA.OELTI.OELTS 
COMMON  /COM«/  FI .W 
COMMON  /COM*/  INO 

EXTERNAL  OIST 

OATA  IERR  /3*(-0).O.t»(-Q)/ 

CALL  SYSTEMC  (3*. IERR) 

CALL  SYSTEMC  U IS. IERR) 

WRITE  (A. 100) 

READ  INFUT 

INITIALIZE  VARIABLES 

5  AO  (5*100)  N.K.&IOMA. OELTI. DELTS 
IF  (K.LT.3  .OR.  K.OT.N1  00  TO  IB 
IF  (N  .(IT.  30)  WRITE  It  .100 )  N 
IF  (OELTI  ,BT .  I0.PSIBMA)  OELTI *1 0.*S IOHA 
IF  (OELTI  .BT,  10.P1IOMA)  DELTIwIO ,»S I OMA 
RCSlwO. 

FESIaO. 

PI ■]. I*  I502AS3A 
NN*N/B*B 
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uuu  UUU  UUU  O  OUO 


I f  (N  .EQ.  SO)  NN«1Z 
ERROR*10.»a<-NNI 

ab«io.*iioma 
Ra*10.  •HOB* 

EVALUATE  TRIPLE  INTEGRAL • 

IN0»l 

PCSl aSQUANK  ( A. 0. ERROR. RUM. OIST) 

If  IN  .EQ.  K)  00  TO  T 

IND*Z 

RESZ*SQUANK  ( A 1 1 .ERROR .RUM.0I*T) 
DETERMINE  FACTORIALS 

T  CONTINUE 
NfACTal 
KZPACTal 
KSfACTal 
NKFACTat 
NKIFACTal 
00  10  I-l.N 
NfACTaNFACTM 


If 

(I 

•  LE. 

(K-Z) I 

KZFACT  aKZFACTal 

If 

<1 

.LE. 

<K-3)  ) 

K3FACTbK3F  ACT* I 

If 

(I 

.LE. 

(N-K) > 

NKFACTbNKFACT*! 

If 

(I 

.LE.' 

(N-K-D)  NKlfACTaNKlf ACT»I 

10  CONTINUE 

COMPUTE  AND  PRITC  PR01A0ILITIES 

PROS 1"NF ACT/ <NKFACT»K3FACT) *RES1 
RflO»|»NFACT/ <NK1FACT»KZFACT) aRESI 
PRQBaPR0f)l*PRO9t 

WRITE  (0.300)  N.K.5IQMA.DELT1 .DELTZ.RR0B1 .PROBZ.PRON 
00  TO  S 

19  aRlTE  (0.4001  N.K 
STOP 

100  FORMAT  (ZI9.3FZ0.7) 

ZOO  FORMAT  (1H1//) 

300  FORMAT  (1H  »4MN  a  .12.3X.4HK  a  . IZ.3X .OHSIOMA  a  .F10.T.3A. 

•  0H0ELT1  a  .F15.0.3X.0HOELTZ  a  • F 19.0 • 3X/1H0 . 

•  3THPR0BA11LI TV  (1  THROUOM  K  FUNCTION)  a  .F10.6/1M  » 

•  39HPROSAOILITY  (Z  THROUOM  K*1  FUNCTION)  a  .P0.6/IH  , 

•  ZZHPRORARILITT  (TOTAL)  a  ,F2 5.6/////) 

*00  FORMAT  (in  . 3SH*a**»  INVALID  VALUE  OF  N  OR  K  •♦•••tJKt 

•  *HN  a  .I2.3X.4HK  a  ,13) 

900  FORMAT  I21H  *»••*  WARNINU  -  N  a  ,JZ. 

•  91M.  CHECK  TmE  TOLERANCE,  LIMIT*  OF  THE  INTEGRALS  •*••*///> 
END 


nnnon  non 


FUNCTION  OUT  (Ml 


C 

c 

c 

c 

c 


c 

c 


c 


c 


COMMON  /CONI/  N«K*S10MAi0ELTI*0CLTS 
COMMON  /COMJ/  MItNM 
COMMON  /COM*/  IND 


external  rxx.rxy 


nN*n/9*6 

If  (N  .CO.  90)  NN-12 
T0LX«10 •••  (**NN) 

TOLY"  l  0  (”NN) 

RHlM"l«/(SORTU."H>MlOMA|«EXR< 


•I. 


•«••{/ (2.*SI9MA**Z) > 


If  (INO  .CO.  *)  00  TO  5 
OOMNnO, 

UMOILTl 
CO  TO  10 
9  OONN"OELT| 

If"*0«tTlTluT.  IO.MIOMA)  0R"U.«SI#MR' 

10  CALL  OOLINT  (UOMNiUR.F XX»FXY«TOLX*TOLY.AN*.RUN.AUHYM) 

oist«i»hi**an* 
return 
I  NO 


NORMAL 


. 


C 


c 

c 

c 


SUBROUTINE  /X*  (V.Vl.ri.VJI 


THU  ROUT  INK  «BOVIOM 
AS  FILL  At  TMC  LIMITS 


ThC  INTtORANO 
Of  INTtORAT ION 


FOR  THE  BCCONO  INTIORAL 
rOR  THE  FIRST  INTIORAL 


COMMON  /COM!/  N,K.SIONA»0ELTl«DELT« 
COMMON  /COM*/  Rt»M  „ 

•HlYMMl  •  /  UORT <t«*RX I *9 IONA) *EXRI"i •  * 


(V*»»*F»/«I.#SIOMA»«l) > 


NORMAL 


Yl*FWlV" 

Y2»0. 

Y3"OELT2 

RETURN 

END 


non  n  n  n  n  o  nnnn 


JUNCTION  PAY  IV.UI 


THIS  ROUTINE  PROVIDES  THE  INTEONAND  POR  THE  FIRST  INTEGRAL. 

COHMON  /CONI/  N.K.SISMA.OELTl. CELTS 
COMMON  /CONI/  RIiM 
COMMON  /COMA/  I NO 

RHIUVUaI ./ (SORT (S.*Rl)*BISMA)*EXPIsl,*(U*V*M)A*|/ (Z,*SIQMA»*|) )  NORMAL 

CAPVWaPNO  < <V*»I /SISMAI  NORMAL 

CAPUVMaPNO  I (U*V*R) /SISMAI  NORMAL 

IP  UNO  .SO,  |)  SO  TO  19 

IP  (N  .EQ.  SI  SO  TO  I 
IP  (N  .EQ.  K)  SO  TO  9 
IP  (K  ,EO,  SI  SO  TO  10 

PXYa tCAPuVV-CAPVR I •• |K»|) *PMIUV<* ( l ,-CAPuVM) •*<N-K I 
SO  TO  IS 
I  FAYaPHIUVM 
SO  TO  ZS 

S  FXValCAPUVV-CAPVW)  ••IK-1)  MMIUVR 
SO  TO  SB 

lo  FXY«PH!UVM*ll.-CAPUVtf)*«<N»K> 

SO  TO  IB 

IB  CONTINUE 

IP  (N  .EQ.  (K*l|)  SO  TO  ZO 

PXVa  ICARUVV-CAPVMl  ••  IK-1)  ARHIUVO 1 1  ,-CAPUVW I  •*  I  N-K-l  I 
SO  TO  ZB 

10  PXV»  ICARUVV-CARVMI •• <K-*| APMIUVK 

IB  RETURN 
END 
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